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Abstract 
(N, 

K^ , In this article, we investigate the Drell-Yan process of the hght neutrahno pair x^x'j {^i3 = 1) 2) 

CN . 

t^ \ productions at proton-proton cohisions and we present the general formulae for the differential 



cross sections. We conduct an extensive examination of the dependence of the total cross section 



o 

Q . of the subprocesses qq — > x?Xj o^i the beam energy, on the mass of the squarks and also on the 

i-^ \ M2 gaugino for the three extremely different scenarios. For all three cases, the outcomes are as 

Q_i' follows: The dependence of the total cross section of the subprocesses qq -^ Xi'x? on the beam 

(D : 

rS^ • energy is dominated by one of the subprocesses, qq -^ XiX2- O^ the other hand, the dependence 

of the total cross section of the subprocesses qq -^ X^ix'j on the mass of the squarks is dominated 
5h ' by one of the subprocesses, qq — > XiXi- We derive there from that our findings may lead to new 

insights relating to experimental investigations and these dependencies may be used as bases of an 
experimental research for the neutralino pair at LHC. 
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I. INTRODUCTION 

The Standard Model (SM) is a successful theory of strong and electroweak 
interactions up to the energies accessible at present [1]. The hierarchy prob- 
lem suggests that, in principle, SM is one of the fundamental effective theories 
of the low energy region. Supersymmetry (SUSY) is presently the most pop- 
ular attempt to solve the hierarchy problem of SM, where the cancelletion of 
quadratic divergences is guaranteed and hence any mass scale is stable under 
radiative corrections. The most favorable candidate for a realistic extension 
of SM is the minimal supersymmetric standard model (MSSM). In MSSM, 
a discrete symmetry called R-parity [2,3] is kept in order to assure baryon 
and lepton number conservations since the gauge-coupling unification sup- 
ports conservation of R-parity. The minimal supersymmetric standard model 
(MSSM) [4] predicts that there exists an absolutely stable LSP. Most often 
the LSP in the MSSM theory is the lightest Majorana fermionic neutralino 
X^. Therefore the production of the lightest neutralino x^ and the second 
lightest neutralino X2 ^^^Y be studied at present and future experiments and 
the detailed study of the neutralino sector will help us to determine which 
kind of the supersymmetric models really exists in nature. They are de- 
termined by diagonalizing the corresponding mass matrix. In MSSM, the 
mass matrix depends on four unknown parameters, namely fi, M2, Afi, and 
tanP = v^jvi^ which is the ratio of the vacuum expectation values of the 
two Higgs fields, ^x is the supersymmetric Higgs-boson mass parameter and 
Ml and M\ are the gaugino mass parameters associated with the SU{2) and 
U{1) subgroups, respectively. The direct search of supersymmetric particles 



in experimental research is one of the promising tasks for present and future 
coUiders. The multi-TeV Large Hadron CoUider (LHC) at CERN and the 
possible future Next Linear CoUider (NLC) are elaborately designed in or- 
der to study the symmetry-breaking mechanism and the new physics beyond 
SM. If the supersymmetry really exists at TeV scale, SUSY particles should 
be discovered and it will be possible to make accurate measurements to de- 
termine their masses and other parameters of the Lagrangian at LHC, and 
then we will have a better understanding of the supersymmetry model. We 
know that there are several mechanisms inducing the production of a neu- 
tralino/chargino pair at hadron colliders. One is through the quark- antiquark 
annihilation, called the Drell-Yan process, and another is via gluon-gluon fu- 
sion. Although the antiquark luminosity in the distribution function of the 
proton is much lower than gluon, the cross sections of the neutralino pair 
productions via the Drell-Yan mechanism are competitive with those from 
the gluon-gluon fusion, since the former mechanism of the neutralino pair 
productions is accessible at the tree level. These facts make the production 
rates in the Drell-Yan process competitive with or even larger than those 
in gluon-gluon fusions. But, the reactions qq -^ xfXj are only subprocesses 
of the parent pp hadron collider. In work [5], the authors have considered 
the production of neutralino pair at a high energy hadron collider, putting 
a special emphasis on the case where one of them is the lightest neutralino 
Xi, possibly constituting the main Dark Matter component. Neutralino pair 
production in proton-proton collisions have been studied in [6] as well, but, 
our results disagree with those conducted in Ref. [6]. In our calculations, 
we have used the anticommuting nature of the Fermionic fields in amplitude 



and cross section of processes, which do not agree with the fermionic symme- 
try property assumed in Ref. [6]. Therefore, this approach is significant for 
the theoretical and experimental studies relating to the neutralino pair pro- 
ductions through the proton-proton collisions at LHC. These results imply 
an interesting complementarity between the future LHC measurements and 
the related 77 -^ X^X^ measurements at a future Linear Collider. Within 
this context, this paper is organized as follows: in section Ull we present 
some formulae for the neutralino/chargino sector. In section Hill we provide 
the formulae for the amplitudes and the differential cross sections of subpro- 
cesses qq -^ XaX^ and in section IIVI we present the numerical results for the 
cross-section and discuss the dependence cross-section on the SUSY model 
parameters. We state our conclusions in section IVl 

II. MSSM PARAMETERS IN NEUTRALINO/CHARGINO SECTOR 

In the minimal supersymmetric extension of the Standard Model (MSSM), 
the physical neutralino mass eigenstates Xi (i=l52,3,4) are the combinations 
of the neutral gauginos ( B and W^) and the neutral higgsinos {Hi, H2). 
In the two-component fermion fields iJj'j = (— ^A-^, —iX^^ipHo^i/jjjo) [4,7],where 
A^ is the bino and A^ is the neutral wino, the neutralino mass term in the 
Lagrangian is given by 

L= --{i/jYMilj^ + h.c. 



The neutralino mass matrix [4,7] in the {B, W, H^, H2) basis, 

'Ml -mzcpsw mzspsw ^ 



M 



M2 mzcpcw -mzspcw 

-mzcpsw mzc/scw -fi 

\ mzspsw -mzspcw -/i ^ I 

is built up by the fundamental supersymmetry parameters: the U{1) and 
SU{2) gaugino masses Mi and M2, the higgsino mass parameter fi, and the 
ratio tanP = 'U2/'yi of the vacuum expectation values of the two neutral 
Higgs fields, which break the electroweak symmetry. Here, SjS = sin(3^ cp = 
cos (3 and s^k, cw are the sine and cosine of the electroweak mixing angle 
9w- In CP-noninvariant theories, the mass parameters are complex. By the 
reparametrization of the fields, M2 can be taken as real and positive without 
loss of generality so that the two remaining nontrivial phases, which are 
reparametrization-invariant, may be attributed to Mi and /i: 

Ml = iMile'-^i and /i = \ii\e''^^ (0 < 0i,0^ < 2%) 

The experimental analysis of neutralino properties in production and decay 
mechanisms will unravel the basic structure of the underlying supersym- 
metries theory. The charginos xt {J — I5 2) mass matrix in the current 
eigenstate basis have the form [7] 



M,= 



M2 V^mwCfi \ 



is diagonalized by two different unitary matrices UuMJJ^ = diag {m^, m^}. 
parametrized in general by two rotation angles and four phases: 

\ -sl cl 
and 

\ -sr Cr 

where cl^r = cos(J)l^r and sl^r = sin(j)L,Re^^^'^ ■ In the limit of M|,|;Up ^ 
m| and IM2 ± {nW^ ^ m|, the following expressions 

mf = M2 + X2[M2 + |Ai|s2/3COS0j, 
^2 = ImI^2[|m| + ^2S2/3COS(/)J, 

are found for the chargino masses and 

^L = T71 my^^P + ^ ^P) ^1 = +A2--pS2/3SZn0^ 

iK/2 — |/i| ^VJ-2 

Sr = V72 ^(/^C/3 + M2S/?) 72 = -X2-r—rS2pSin(f)^ 



2 2 



|M2|2-|^|2 

for the mixing angles and phases. In the present work, we have investi- 



gated the Drell-Yan process of the hght neutralino pair XiX] {hj = l?^) 



J 



productions at hadron colliders. Since the neutralino mass matrix M is sym- 
metric, one unitary matrix A^ is sufficient to rotate the gauge eigenstate basis 
(B^, W^, Hi, H2) to the mass eigenstate basis of the Majorana fields xl*- 



Md = N^MN = ^m^uEj, (2.1) 



where A^ is a unitary matrix. To determine N, it is easiest to square 
eq.(2.1) obtaining 

4 

Ml) = N-^M+MN = ^ m|o^j, (2.2) 

where {Ej)4x4 are the basic matrices defined by {Ej)ik = ^ji^jk and x? stand 
for the four component Majorana neutrahnos: 



x? = 



^x°^ 



\ x° / 



,j = l,....,4 (2.3) 



Here, we suppose that the real eigenvalues of Md are ordered in the following 
way 

Xl A2 A3 Xi 

The mass eigenvalues m^o (j = 1,2,3,4) in Mjj can be chosen as positive by 
a suitable definition of the unitary matrix A^. In this work, we consider the 
higgsino/gaugino sector with the following assumptions: First, for simplifi- 
cation, CP-conservation is hold, namely (f)^ = (f)i = 0. The physical signs 
among Mi, M2 and fi are relative, which can be absorbed into phases 0^ 
and 01 by redefinition of fields. Thus, Mi, M2 and fi are chosen to be real 
and positive, i.e., Mi, M2, fi > 0. With the above assumptions, there are 
several scenario for the choice of the SUSY parameters for the investigation 
of the neutralino pair production in hadron collider. One can employ the 
scenario of taking Mi, M2, fi, and tan/3 as input parameters, and then get 
all the physical chargino and neutralino masses and the matrix elements of 
Ur, U^ and A^ as outputs. Also, there are other alternative scenarios, such 



as the CP conserving mSUGRA scenario with five input parameters, namely 
^1/25 ^0, ^0, A^ and tan/3, where mi/2, tuq and ^40 are the universal gaugino 
mass, scalar mass at GUT scale and the trilinear soft breaking parameter in 
the superpotential respectively. From these five parameters, all the masses 
and couplings of the model are determined by the evolution from the GUT 
scale down to the low electroweak scale [8]. Since SUSY parameters should 
be extracted from the physical quantities, one can also choose an alternative 
way to diagonalize the mass matrix M, by taking any two physical chargino 
masses together with tan/3 as inputs. There are several scenarios about the 
choice of two chargino masses and tanP [10]. Also there are two possible 
scenario about the choice of tanP: scenario with small tan(3{tan(3 ^^ 1 ^ 3) 
and scenario with large tan(3{tan(3 « 30 -^ 70) [9]. In this work, we take 
two chargino masses m~+ and scenario with small tan(3 as inputs. In this 
way, the two fundamental SUSY parameters, M2 and ^ can be figured out 
from the chargino masses by using the following formula: For given tanf], 
the fundamental SUSY parameters M2 and fi can be derived from these two 
chargino masses [11]. The sum and differences of the chargino masses lead 
to the following equations involving M2 and fi: 



Mi + \ii\' = m~+ + m^+ - 2miy, (2.4) 



Mi\ii\^ - 2m^sin2/3cos0^M2|^| + (m^sm^2/3 - mi+mi+) = 0. (2.5) 



Xi X2 



The solution of (2.5) is given as: 

^2|a^| = 'm-^cos(l)^sin2(3 ^ ^/m1+m'i+ — m'^sin'^2f3sin'^(j)^. (2-6) 

From (2.4) and (2.6), one obtains the foUowing solutions for M2 and fi: 

2M| = (m|^ + m|+ - 2ml^) T (^(m|+ + m|+ - 2m2^)2 - A±) , (2.7) 



2|^|2 = (m|+ + m|+ - 2ml) ± (^(m|+ + m|+ - 2m2J2 - A±) , (2, 
with 



A± = 4 



m~+m~+ + m^^cos2(j)^sin 2(3 ± 2my^cos(j)^sin2f3- 



> 



^/m1+m1+ — rni^si'n?2(3si'n?c 
Y Xi X2 

where the upper signs correspond to M2 < |/u| regime, and the lower ones to 

M2 > \fi\. Therefore, for given tanP, M2 and fi can be determined in terms 

of the masses of the charginos m~+ and m~+ by using (2.7), and (2.8) from 

which one gets four solutions corresponding to different physical scenarios. 

For fi < M2, the lightest chargino has a stronger higgsino-like component and 

therefore is referred to as higgsino-like. The solution fi > M2, corresponding 

to the gaugino-like situation, can be readily obtained by the substitutions: 

M2 -^ fi, and fi -^ sign{fi)M2. In this paper, we assume the GUT relation 

11,12] 

5 
Ml = -M2taT?Bw- (2.9) 



Thus, the neutrahno masses can be determined by solving the characteristic 
equation associated to this system, that is 

X^ - aX^ + 6X2 - cX + d = 0, ^2.10) 

where 

a = Ml + 2/1^ + Ml + 2m|, 

xcos(l)ij, — 2m%s'^Misin2f3cos{(j)^ + 0i), 

c = 11^ Mf + ii^m\sin^2(5 + Mfm|c^(2;u2 + m|c^)+ 

^l(^|svF+V(^|s?^+^?)+A^^)-2/im|s^Mi(^2^M|)sm2/?cos(0^+0i)+ 

2m|c^M"2[m|Mis^cos0i — ii{ii^ + Ml)cos(l)^sin2l3], 
d = m^Cy^rfi M^sin 2(3 + 2111^^1 MiM2Cy^r(mzSy^rsin2(3cos(j)i — fiMicos(j)^)+ 

li^m\s\rMlsin2(5{m\s^sin2(5 - 2iiMicos{(f)i + 0^)) + j/mImI. 
Solving Eq.(2.10), we get the exact analytic formulae for the neutralino 



masses 



2 2 a / _ 1 I ¥ 



2 2 a , /_ 1 /I , V 



Kl'Kl = i+2^2V'"'' + 47' ^^-^^^ 



where 



"^ ~ (3 • 21/3) + 3 • g 
Q = (A;+ v'A;2-4/i3)i/3^ 

10 



k = 2b^ - 9abc + 27c^ + 27a^d - 72bd, (2.12) 

h = b'^ -3ac+ 12d, 
p = a^ — 4ab + 8c, 

Starting from Eq.(2.2), we get 

{M+M)N - NMl = 0. (2.13) 

A more explicit form of this matrix equation is 

(^11 - m|o)iVi, + AuN2j + Ai^Nsj + AiMj = 0, 
A21N1J + {A22 - mlo)N2j + A23NSJ + A24N^j = 0, 

AsiNij + ^32A^2j + (^33 - mlo)Nsj + ^34A^4, = 0, 

A^lNij + ^42A^2, + ^43A^3, + (^44 " m|o)iV4, = 0, (2.14) 



11 



j=l,...,4, where A,, = ELi ^^W^r 

All = Ml + m|s^, 
Ai2 = A21 = -mlswcw, 
Ai3 = A31 = -MimzCfjSw - iimzswsp, 
Au = A41 = Mimzspsw + /imzCissw, 
A22 = Mi + m|c^, 
-423 = -432 = M2mzcpcw + jimzspcw, 

A 2,22 

Aii = ji + ruzCp, 

^24 = -442 = -M2mzspcw - idmzcpcw, 

Am = -443 = -m%spci3, 

A 2 2,2 

^444 = ^^S/3 + A^ • 

The diagonahzing matrix A^ can be obtained by computing the eigenvectors 
corresponding to the eigenvalues given in Eq.(2.11). Indeed, by inserting a 
generic eigenvalue m^o, into Eq.(2.14) and dividing each one of these equa- 
tions by A^ij, where it is assumed that Nij ^ 0, we get 

A ^^^ ^A ^''^ ^A ^'^ m^ - A 

{A22 - m|o -^ + A2,-^ + A2,-^ = -A21, 
^3 Nij Nij Nij 

Solving this system of equations, and taking into account the relation 

l^^ijf + l^^2/ + \Nsj\' + |iV4jf = 1, (2.16) 

12 



it yields the Nij matrix's component 



lA 



ijl 



<ij 



^^J Vl^bf +l^2,f +|A3,f +|A4,f' 



(2.17) 



when, ^ = 1, ...4. Here, 



Ai, = 



^22 - mlo A23 

^3 



A 



24 



^32 
^42 



-433 - w|o -434 



A 



43 



^44 - m|o 



2, 3, 4, is formed from Aij by substituting the (i-l)th column by 



and Aij, i 
( \ 

-A21 
-Ail 

\ -^41 / 
III. CALCULATION OF THE CROSS SECTION 

The Neutralino pair productions, which can be produced via the collisions 
of quark and antiquarks in protons, can be expressed as 



q{pi)q{p2) -^ x?(^i)x?(^2), 



(3.1) 



where pi and p2 represent the momenta of the incoming quark and antiquark, 
and ki and /c2 denote the momenta of the two final state neutralinos, respec- 
tively. The Mandelstam invariant variables for subprocess (3.1) are defined 
as 

s= {pi^ P2f, i = {pi - kif, u= {pi- k2f. (3.2) 

The Feynman diagrams of the subprocess are shown in Fig.l. The relevant 
couplings of the supersymmetric particles are deduced from the following 



13 



interaction Lagrangians of the Supersymmetric Standard Model [13]: 

Lz^^^ = \-^Z,tl'{OU,PL + 0:;.i?,P^)x?, (3.3) 



Lz^qq = —^qi^{L,PL + RqPR)qZ^, (3.4) 



Lqq^. = -V2gq[af{qn)PL + af (gn)i^i?]x?9n. (3.5) 

In Eqs. (3.3-3.5) Xa ^^^ Q ^^^ four-component spinor fields and q is the field 
of the squark. Furthermore, g = e/sinOw (e > 0) is the weak coupling 
constant, Qw the Weinberg angle, Pr^l = |(1 ± 7^), while the coupling 
constant Of^-, L^, Rq and a- ' {qn) are given by 

0:f = \{N,,N;, - N,,N;,)cos2f3 - hN,,N;, + N,,N;,)sin2/3, (3.6) 



o;f = -o:f , (3.7) 



Lq = 2/3(1 - 2szn2eT4^|Q,|), Rq = -2sin^ewQq. (3.8) 

with /^, Q„ being the isospin and charge of the various g^-quarks, and 



'q^ ^q 



afiuL) = -— ^= {Nusw + SN2iCw), 



3y'2swCw 



e 



<(^i?) = 1^ K. 

oCw 



afidi) = - r- {'Nusw - 3iV2zCvi/), 

6\J2swCw 

6Cw 
14 



T /^"^ \ ^''^U TV 7- 

a, [Ur) = ^= N4^, 

ylmwswsp 

\/ZmwSwS(3 
afidn) = - - """ N,., 



af(d,) = - ^ "'"■' iV|, (3.9) 

In (3.9), (g = u,d) refer to the incoming up and down quark (antiquark) 
of any family, while {qn = QL^qR) denote the corresponding squarks. We 
also note that the mixing matrices N in (3.6, 3.9), control the Bino, Wino, 
Higgsino components of the neutralino in the ^x^x^ and qq)('^ coupling. The 
corresponding Lorentz invariant matrix element for each of the diagrams can 
be written as 

T = r, + r£ + r^, (3.10) 

where 

7) = ^ 7 —Ui{ki){af{qn)PL + a^{qn)PR)u{pi)v{p2)x 



^ 9 

t — m~ 

Qn 



{aY{qn)PL + af{qn)PR)Vj{k2), (3.11) 

Tu = -Y^ ^ _ ^ Uj{k2){af*{qn)PR + af*{qn)PL)u{pi)v{p2)x 

Hi If t^- 

n Qn 

iaf{qn)PL + af{qn)PR)v,{ki), 

where the index n refers to the summation over the exchanged L- and R- 
squarks of the same flavor in the t-and u- channel, and (z, j) describe the final 

15 



neutralinos. From the total (gauge-invariant) amplitude T, which is the sum 
of the partial amplitudes Eq.(3.11), we obtain the differential cross section 
as 



'^'' - ^'' ' "^ ^^" {Mrs + Mii + M^^ - 2M,i + 2M,^ - 2Mi^) , (3. 12) 



dQ. 3847r2s2^2 
where 



Kj = ./(s - m|o - m|o)2 - 4m|om|o/2, (3.13) 

y A 2 Aj' Ai A J 

and {\Y'^ is the final identical-particle factor. The squares of the matrix 



-2 

element have the form as 

A 



4 sm^ 6^ cos^ 6vK y y x, x^ 



mlo - i){m\o - i) - m^om^os(Of + Of\ (3.14) 

-Vi Aj 

1 



^ii = 7? T^7 TT(af (^fc)af*(^/) + af{qk)af\qi)){a^{qk)aY{qi) + 



(^-%J(^-%) 



af{qk)af\qi)){m\. - f)(m|o - f), (3.15) 



' i*/~\„i/~\ I „i?*/~ \„i?/~\\/„i/~\„i*/ 



M^^ = t:^ 2^7:^ 2\(^i{qk)at{qi)+af*{qk)af{qi)){a^{qi)aj*{qk) + 



afiqi)af"iqk))imlo - ii)(m|o - u), (3.16) 



-J v-ity-^j K^n^^v^Xt "^^v"% 



^^t^ = 7T YTT^ -Al [af"{Qk)a^{Qi)af{qk)a^*{qi)+a^*{qk)af{qi 



16 



af''{qi)a^{qk)] ((mL-u)(rn|o-w)+(m|o-t)(rn|o-f)-s(s-m|o-m|o))+m|om|o 

•J -■ A 7 A.i A.J Ai Ai At Af A 7' 

xs[af (gO«f (^fc)«f (^/)«f (^fc) + af (g/)af (^fc)af (^0«f (^fc)]} (3-17) 

M,u = ,^-^ ^^^j-. -\{Re[Dzm)[L,af^{qk)aj{qk)Of - 

Rqaf\qk)af{qk)0'^]{m\^ - w)(m|o - u) + [Rqaf\qk)af{qk)Of - 



Lqaf*{qk)a^{qk)0'^]m^om^os L (3.18) 



e 



,2 
Zsin'^Ky]^cos'^Ky-^[t — m^ ^ ' 



Mu = --^ ^7^-^ ^ (i?e[Dz(s)])[i?,af (^fc)af(^fc)0; 



^*^;r,A^^^;r, Ar»*i* 



gfc^ 



L,af (gfc)af (^fc)Og](m|o - f)(m|o - f) + [L,af (gfc)af (gfc)Of - 

i?,af(gfc)af(gfc)0|]m^om3^os|, (3.19) 

The following abbreviation has been used 

Dz{s') = ,\. ^ ■ 

s^ — m^ + imzi- z 

For calculation, we assume uiz = 91.1887 GeV and the widths of the gauge 
boson by Tz = 2.499947Gey. The basic parton model expression for the 
hadron-hadron collision hi{pi)h2{p2) -^ )(A{ki)x^j{kj), [14,15] is 

da{hi{pi)h2{p2) -^ xl{ki)X^j{kj)) = ^ dxidx2Gq^/h,{xi,Q)Gq^/h^{x2,Q)- 



daiqm ^ X?X?)t-^ (3.20) 

-L 1 Oqiq2 

with x^, Xj being the two produced massive particles of mass m^o, m^o. Here 
Gqjhiixi, Q) is the distribution function of partons of type {qi = q,q), in the 

17 



hadron of type hi at a factorization scale Q. Taking the hih2-c.m.. system as 



the lab-system, the lab-momenta of the produced x^ and x^ are [16] 

/cf = {Ei, kr, kiCosO), k^ = (E'j, — /cr, kjCosO), 
where their transverse momenta are obviously just opposite 



(3.21) 



kx — kx. — —kr,, 



(3.22) 



while their transverse energies Et^ 
to define 



Jk^ + mlo, Et^ = Jk"^ + m|o are used 



2E' 



XT,. = 



T,, 



Pt, = kr/Er, = 



\ 



4m~o 

Xi_ 

2 ' 



sx 



T, 



(3.23) 



Note that 



2E' 



XT, 



T, 



V~s 



, A- = kr/ET, = 



\ 



4m~f 



sx 



2 ' 



(m|o - m|o) 

rf _ 172 , _2 _2 _2 _ _2 , /I >^j X» 



"T, 



^T, + "^y" - "^y". ^T, = Xt, + 4 



(3.24) 



(3.25) 



The rapidites and production angles of x^, Xj? ^^ the lab-system, are related 
to their energies and momenta along the beam-axis of hadron /ii, by 



1 , Ei + kiCosOi 1 , Ej + kjCosOj 

-In— : —, Vj = -In-r -^ r^ 

' ^^ 2 Ej-kjCosOj 



The center-of-mass rapidity y of the Xi'Xj pair, and their respective rapidities 
y* in their own cm. frame, are defined as 



yi = y + yh yj = y + y1, 



(3.27) 



A = 2/^ - yj 
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yz -y 



J' 



(3.28) 



The fractional momenta of the incoming partons are expressed in the terms 
of their lab-momenta by 

1^1 = 9(^1'^'^' ^l)' 1^2 = 2(^2,0,0,-0:2), l? = l?l +j32, (3.29) 



fs fs 

V^ = —(3^1 + 3^2) = E^^Ej, ps = — (X1-X2) = {hcos9i + kjCos9j), (3.30) 

which lead to 

XI = lixT^ey^ + XT,e^^] = ^e\ (3.31) 

2 Vs 



X2 = l[xT^e-y^ + XT^e^y^] = ^e-\ (3.32) 



s = M^ = {pi + j92)^ = 3:1X25 = -[xy^ + Xy^. + 2a:r,a:T^cos/i(A?/)]. (3.33) 

Using this, s, xi, 0:2 may be calculated in terms of the final particle rapidities 
yi, Hj and their transverse momenta. From them, y is also obtained through 
(?/*, ?/*). The remaining Mandelstam invariants of the subprocesses satisfy 

i = (j9i - kif = m|o - M(E^ - k^cosO*) = m|o - ^Mv^e"^^' = (3.34) 

"^|o - ^xiXT.e-y- = m|o - M(E* - k^cosO^) = 
At 2 ^3 •' 



m~o 



2 ^T, ^^J /-y* _2 ^ 



Mv^e^^* = m|o - -X2Xr,e^% (3.35) 



M = (j9i - A; .)2 = ^2 _ ^^^* ^ rcos^*) = m|o - ^M^ie^^* 

A 2 Ai ) 

m| X2XTe^' = m|o - M(E* + k*cos6*) = 
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m~o 



2 '^Tj n ^ r- —II* 9 S 



M^se-y''^ = mlo - -xiXT-e-y^. (3.36) 



^j 2 ^i 2 



r = - = xia;2, (3.37) 

s 

where ^* describes x^ production angle in x^x^ -cm. frame (the x^ one being 
TT — ^*). The energies of the two final particles in their cm. -frame are 

- , 2 2 ^ , 2 2 

s + m~o — m~o s + m~o — m~o 
E; = ^V ^, £:* = ^ ^, (3.38) 

their momentum is 

^' = 2lf f^^' " ""^^ " ""1°^' " 4m|om|o]0-5, (3.39) 

and their velocities 

[(M2 - m|o - m|o)2 - 4m|om|o]i/2 

A* = r/E* = j^^-T^ W^^^ (3-40) 

[(M2 - m|o - m|o)2 - 4m|om|o]^/^ 

/?; = r/E; = J; "- "' "- , (3.4i) 



We also have 



co.r^'""^^*- ^-^. 



(^t /?; 



smr = ^, (3.42) 



^.-,^yt-^_Z^-l±P>^ (3 43) 



u — m~ 1 + 6*cosO* 

Xj J 
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/3>sr = ^ = ^, (3.45) 



Xz(m|. -m|) +M2 
4. - "^1 



X^- = ^^'\/2\ 2^^ -^ (3-46) 



(M2 + m|o - m|o)'x* - M^mloil + X.)' 

4 = ,1.^ ^2 ' (3-48) 

M^(1 + Xj)^ 

4(M2 + m2 - m 0)2^, 

4- - ^^""^tTf-'f"-. (3.50) 

^^ M2s(l + xj)^ 

Using (3.20) we define tlie expression for tlie cross section in terms of the 
overall center-of-mass rapidities of the two jets yielding 

-J = XiX2^Gq^/h,{xi, Q)Gq^/h^{x2, Q)-j{qiq2 ^ X?Xj)- (3.51) 



dy.dyjda ' ^^ '''"'' ''-' '''"'' ^'-'di 



IV. NUMERICAL RESULTS AND DISCUSSION 



In this section, we discuss the numerical results for the process pp -^ X^X? 
at LHC (cm energy -^5=14 TeV) generated by the subprocesses qq. As we 
assume, Xi is likely to be the LSP, the three types of channels: qq -^ XiX?, 
qq ~^ X2X25 QQ ~^ X1X25 would be the most dominant neutralino pair pro- 
duction processes, which may lead to the first detection of SUSY particles at 
the LHC. The numerical results of the cross sections of these three processes 
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have been presented and their dependencies on the basic SUSY parameters 
have been discussed. We divide the input MSSM parameters into two parts. 
One part is for the general parameters included also in SM, and the other 
part is the ino and squark sectors of MSSM. For the first parameter part, 
we take m^o= 91.1887 GeV, sin 6^ = 0.2315, a^=l/137. For the second 
part, we just limit the values of Mi, M2 and fi to be real, positive and below 
1 TeV, and take tan/3 = '2,mui — ^^=300 GeV, rriu^ = m^=500 GeV. 
Also, we fix the heavy chargino mass as m~+ = 4500 eV and for the lightest 
chargino mass m~+ = IbOGeV. By using Eqs.(2.7, 2.8) with above chargino 
mass values, one may have two choices of parameter sets for /i and M2 in 
two extreme cases, which are the Higgsino-like and the gaugino-like respec- 
tively. For the Higgsino-like case, we get M2 = 437 .96GeV , n = 169.7b3GeV, 
Ml = 219.703Gey and by inserting the values of M2, /i and Mi into Eq.(2.11) 
for neutralino masses, we get 

m^o = 123.242GeV, m^o = 179A7bGeV, m^o = 238.199Gey, m^o = 4b8A79GeV 

For the gaugino-hke case, we have M2 = 169.753Gey, ^ = 4:37.96GeV, 
Ml = 8b.lb7GeV and by inserting the values of M2, /U and Mi into Eq.(2.11) 
for neutralino masses, we then get 

m^o = 77.212Gey, m^o = 153.859Gey, m^o = U9A20GeV, m^o = 4b2A43GeV 

For full discussion, we present the results for the mixture case as well. We 
take the heavy chargino mass m~+ = 280. GGe^ and for the lightest chargino 
mass m~+ = 128GeV, the corresponding outputs obtained as M2 = fi = 
200.598Gey, Mi = 100.63Gey and also by inserting the values of M2, ^ and 
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Ml into Eq.(2.11) for the neutralino masses, we get 

m^o = 72.b76GeV, m^o = U7.127GeV, m^o = 207.184Ge"\/, m^o = 278.63bGeV 

Xl ' A2 ' A3 ' A4 

As an example for the quark distribution function inside the proton, we use 
the MRST2003c package[17]. After this, the cross sections of the subpro- 
cesses qq -^ XiX?, X1X2' X2X2 can be numerically evaluated. For illustration, 
we have calculated the total cross sections of the subprocesses uu -^ XiX?, 
X1X2. X2X2. and dd -^ XiXi, Xixt X2X2. the dependence on the beam en- 
ergies, on the mass of squarks and also of the M2 gaugino mass. We will 
illustrate this for three extremely different scenarios. In Fig. 2-4, we show 
the dependence of the total cross sections for the subprocesses uu -^ XiX?, 
X1X25 ^X2 c>f the beam energy. In our calculations, the beam energy for the 
subprocesses s changes in the region (400 -^ 5000) GeF, which corresponds to 
the beam energy for process s (500 -^ 16000) GeF. As shown in Fig. 2-4, all 
subprocesses in the gaugino-Hke scenario, the total cross section is 30 per- 
cent larger than the mixing scenario in magnitude, and larger than in the 
Higgsino-like scenario as 1.79-2 order of magnitude. As seen from Figs. 2-4, in 
all three scenarios, the dependencies of the cross section on the beam energy 
demonstrate the same behavior. Also, as seen from Figs. 2-4, all three sce- 
narios, by increasing the beam energy from 400 GeV to 3600 GeV, the total 
cross section is monotonically increasing. In Figs. 5-7, we show the depen- 
dence of the total cross sections for the subprocesses dd ^ XiXv X1X2' X2X2 
of the beam energy. As shown in Fig. 5-6, all subprocesses in the gaugino-like 
scenario, the total cross section is 8 percent larger than the mixing scenario 
in magnitude, and larger than in the Higgsino like scenario as 2-2.137 order 
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of magnitude. As shown in Fig. 7, in the gaugino-like scenario, the total cross 
section is almost 2 percent larger than the mixing scenario in magnitude, and 
larger than in the Higgsino-like scenario as two order of magnitude. As seen 
from Figs. 5-7, in all three scenarios, the dependence of the cross section on 
the beam energy demonstrates the same behavior. With the increase of the 
beam energy from 400 GeV to 4000 GeV, the total cross section is monotoni- 
cally increasing. The total cross section of the subprocesses uu -^ X?X2 in the 
gaugino-like scenario, appears in the range of 2 to 140 fb and should be ob- 
servable at LHC. It should be noted that, one of the subprocesses, qq -^ x\x\ 
dominates all three scenarios. According to our opinion, it may be used as 
a probe for an experimental search on the neutralino pair. Figs. 8-13 show 
the dependence of the total cross sections for the subprocesses uu -^ XiX?, 
X1X25 X2X2) and dd -^ XiXi, X1X2. X2X2. of the M2. As seen from Figs. 8- 
10, the cross section is decreasing when the M2 gaugino mass is increasing. 
But, as seen from Fig. 11, the dependence of the total cross sections for the 
subprocesses dd -^ XiX?, X1X25 X2X2 ^^ ^^e M2 have another character. As 
shown in Figs. 12-13, the cross section is decreasing with increasing M2 in 
the range of 300 to 400 and it has a minimum approximately at one point 
M2 = 400Gey. After this point, the cross section increases with M2 gaug- 
ino mass increasing. Behavior of the cross section is decreasing when the 
M2 gaugino mass is increasing as expected. Because, we have obtained the 
following relation for the dependence cross section depending on the beam 
energy: a{Higgsino — like) < a{Mixture — case) < a{Gaugino — like). 
This relation corresponds to M2{Gaugino — like) < M2{Mixture — case) < 
M2{Higgsino — like)^ respectively. Therefore, the cross section is decreasing 
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when the M2 gaugino mass is increasing. Figs. 14- 19 show the dependence 
of the total cross sections for the subprocesses uu -^ XiX?, X1X25 X2X2' ^^*i 
dd -^ XiXi, X1X25 X2X2' of ^^6 squarks mass. As shown in Figs. 14-19, in aU 
subprocesses and scenarios, the cross section monotonicaUy decreases with 
an increase in the squark masses. Therefore, an experimental search for the 
neutralino pair at lower values of the squark mass is preferable. In Figs. 20- 
22, we show the dependence of the differential cross sections for the process 
PP ~^ Xi'Xj Q-s a function of the kx transverse momentum of the neutralino 
pair at rapidity yi = yj = 0. As seen from Figs. 20-22, the differential cross 
sections decrease monotonicaUy pursuant to an increase in the kx transverse 
momentum of the neutralino pair. It should be noted that, one of the pro- 
cess pp -^ X1X2 dominates all three scenarios. As mentioned in above, that 
our results disagree with the results in [6]. So, the expression for the cross 
sections in Ref.[6] is wrong. Therefore, our numerical results also disagree 
with the numerical results in [6]. We compared our results with [18]. The 
pair production rate of X1X2 ^^^ gluon-gluon fusion is about few ten percent 
of that via quark- antiquark annihilation at the LHC. The cross section of 
X2X2 vis- gluon-gluon fusion is about few percent of that via quark-antiquark 
annihilation at the LHC. 

V. CONCLUSION 

In this paper, we have considered the neutralino pair production processes 
in proton-proton collisions at LHC (cm energy ^/s = lATeV). In the descrip- 
tion, we have taken into account the subprocess q{pi)q{p2) -^ Xi'(^i)Xj(^2)- 
We have given detail illustrations for the center-of-mass energy, squarks and 
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M2 gaugino masses for three extremely different scenarios. For illustration, 
we have calculated the total cross sections of the subprocesses qq -^ XiXi, 
XiX2' X2X2' ^^cl the dependence on the beam energies, on the mass of squarks 
and also of the M2 gaugino mass. We have illustrated this for three extremely 
different scenarios. Fig. 2-7 show the dependence of the total cross sections 
for the subprocesses qq -^ XiXi, X1X25 X2X2 of the beam energy. As seen from 
Figs. 2-7, in all three scenarios, the dependencies of the cross section on the 
beam energy demonstrate the same behavior. The total cross section of the 
subprocesses uu -^ X1X2 in the gaugino-like scenario, appears in the range of 
2 to 140 fb and should be observable at LHC. It should be noted that, one 
of the subprocesses, uu -^ X1X2 dominates all three scenarios. According to 
our opinion, it may be used as a probe for an experimental search on the 
neutralino pair. Figs. 8-13 show the dependence of the total cross sections 
for the subprocesses qq -^ XiX?, X1X25 X2X2' of the M2. Figs. 14-19 show the 
dependence of the total cross sections for the subprocesses qq -^ XiX?, X1X25 
X2X2 of the squarks mass. As shown in Figs. 14-19, in all subprocesses and sce- 
narios, the cross section monotonically decreases with an increase in squark 
masses. It should be underlined that for all three scenarios, the dependence 
of the total cross section of the subprocesses qq -^ XaX^j on the beam energy 
is dominated by one of the subprocesses, qq -^ X1X2 and the dependence of 
the total cross section of the subprocesses qq -^ X^X^ on the mass of squarks 
is dominated by one of the subprocesses, qq -^ XiX?- These findings may be 
used as a probe in an experimental search for the neutralino pair at LHC. In 
Figs. 20-22, we have shown the dependence of the differential cross sections 
for the process pp -^ XaX^j as a function of the kr transverse momentum of 



^i /\.j 
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the neutralino pair at rapidity yi = yj = 0. In this case, one of the process 
PP ~^ X1X2 dominate all three scenarios. These features make the neutralino 
pair production processes rather interesting for testing the SUSY dynamics 
at LHC. The reason is that they provide tests which will be complementary 
to those addressing the cascade decays of initially produced colored SUSY 
particles to eventually x?, which is here assumed to be the LSP; e.g. studies 
of mass spectra and decay branching ratios [19]. In particular, consistency 
checks should thus become available, allowing the strengthening of possible 
constraints on the validity of specific models. Moreover, in such neutralino 
pair production, the role of the Majorana nature of the final state particles is 
more prominent than in decays involving just one neutralino at a time. Since 
no such states have been observed in the past, it would be interesting to have 
eventually some experimental support of our understanding of the Majorana 
nature. These results imply an interesting complementarity between the fu- 
ture LHC measurements, the related 77 -^ X^X^ measurements at a future 
Linear Collider and the Dark Matter searches in cosmic experiments. 
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FIG. 1: Feynman diagrams for qq -^ X^X^ process. 
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FIG. 2: The cross sections of the subprocess uu -^ XiX? &s a function of Vs. The curves correspond 
to: sohd-Higgsino-hke, dashed-gaugino-hke and dotted-mixture cases, respectively. 
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FIG. 3: The cross sections of the subprocess uu -^ X1X2 ^^ ^ function of Vs. The curves correspond 
to: sohd-Higgsino-hke, dashed-gaugino-hke and dotted-mixture cases, respectively. 
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FIG. 4: The cross sections of the subprocess uu -^ X2X2 ^^ ^ function of Vs. The curves correspond 
to: sohd-Higgsino-hke, dashed-gaugino-hke and dotted-mixture cases, respectively. 
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FIG. 5: The cross sections of the subprocess dd -^ XiX? &s a function of Vs. The curves correspond 
to: sohd-Higgsino-hke, dashed-gaugino-hke and dotted-mixture cases, respectively. 
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FIG. 6: The cross sections of the subprocess dd -^ X1X2 ^^ ^ function of yl- The curves correspond 
to: solid-Higgsino-hke, dashed-gaugino-hke and dotted-mixture cases, respectively. 
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FIG. 7: The cross sections of the subprocess dd -^ X2X2 ^s a function of \/s. The curves correspond 
to: sohd-Higgsino-hke, dashed-gaugino-hke and dotted-mixture cases, respectively. 
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FIG. 8: The cross sections of the subprocess uu -^ XiX?; as a function of M2 with ^ = 450 GeV 
and VI = 1.5 TeV 
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FIG. 9: The cross sections of the subprocess uu -^ XiX2 5 ^^ ^ function of M2 with fi = 450 GeV 
and ^ = 1.5 TeV 
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FIG. 10: The cross sections of the subprocess uu -^ X2X2' ^^ ^ function of M2 with /i = 450 GeV 
and a/I = 1.5 TeV 
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FIG. 11: The cross sections of the subprocess dd -^ XiX?) as a function of M2 with fi = 450 GeV 
and ^/I = 1.5 TeV 
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FIG. 12: The cross sections of the subprocess dd -^ X1X2' ^^ ^ function of M2 with fi = 450 GeV 
and ^ = 1.5 TeV 
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FIG. 13: The cross sections of the subprocess dd -^ X2X2' ^^ ^ function of M2 with fi = 450 GeV 
and \/I = 1.5 TeV 
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FIG. 14: The cross sections of the subprocess uu -^ XiX?) ^.s a function of the squark mass at 
beam energy \/s = 1.5 TeV The curves correspond to: sohd-Higgsino-hke, dashed-gaugino-hke and 
dotted-mixture cases, respectively. 
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FIG. 15: The cross sections of the subprocess uu -^ X1X2' ^^ ^ function of the squark mass at 
beam energy yl = 1.5 TeV The curves correspond to: sohd-Higgsino-like, dashed-gaugino-hke, 
and dotted-mixture cases, respectively. 
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FIG. 16: The cross sections of the subprocess uu -^ X2X2' ^^ ^ function of the squark mass at 
beam energy \/s = 1.5 TeV The curves correspond to: sohd-Higgsino-Hke, dashed-gaugino-Hke, 
and dotted-mixture cases, respectively. 
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FIG. 17: The cross sections of the subprocess dd -^ X^Xv ^^ ^ function of the squark mass at 
beam energy \/s = 1.5 TeV The curves correspond to: sohd-Higgsino-hke, dashed-gaugino-hke and 
dotted-mixture cases, respectively. 
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FIG. 18: The cross sections of the subprocess dd -^ X1X2) ^^ ^ function of the squark mass at 
beam energy \/s = 1.5 TeV The curves correspond to: sohd-Higgsino-hke, dashed-gaugino-hke, 
and dotted-mixture cases respectively. 
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FIG. 19: The cross sections of the subprocess dd -^ X2X2) ^^ ^ function of the squark mass at 
beam energy \/s = 1.5 TeV. The curves correspond to: sohd-Higgsino-hke, dashed-gaugino-Uke 
and dotted-mixture cases, respectively. 
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FIG. 20: The differential cross sections of the process pp -^ XiXi ^^ & function of the kx transverse 
momentum of the neutrahno pair. The curves correspond to: soUd-Higgsino-hke, dashed-gaugino- 
Uke and dotted-mixture cases, respectively. 



10" 



o 



10" 



=^ 



5 10" 

D 




200 400 600 800 

k., , (GeV/c) 



1000 



FIG. 21: The differential cross sections of the process pp -^ X1X2 ^^ ^ function of the kx transverse 
momentum of the neutralino pair. The curves correspond to: solid-Higgsino-like, dashed-gaugino- 
like and dotted-mixture cases, respectively. 
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FIG. 22: The differential cross sections of the process pp -^ X2X2 ^^ ^ function of the kx transverse 
momentum of the neutrahno pair. The curves correspond to: soUd-Higgsino-hke, dashed-gaugino- 
Uke and dotted-mixture cases, respectively. 
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